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1.  INTRODUCTION 


The  objective  of  the  contract  "Monte  Carlo  Simulation  of  Ion  Collection 
by  a Rocket-Borne  Mass  Spectrometer"  has  been  to  establish  a connection 
between  the  properties  of  the  ambient  ionosphere  and  the  strength  and  chemical 
state  of  the  flux  entering  a rocket  borne  mass  spectrometer.  The  primary 
concern  for  this  study  was  the  description  of  the  collection  of  positive  ions 
by  a flat  faced  cylinder,  approximately  15  cm  in  diameter,  traveling  through 
a weakly  ionized  gas.  Details  of  the  Monte  Carlo  direct  simulation  procedure 
to  predict  the  collection  of  positive  ions  were  described  by  Vogenitz  (1973) 
and  Sugimura  and  Vogenitz  (1973,1975)  and  will  not  be  repeated  here.  The 
results  presented  include  a detailed  examination  of  positive  ion  collection 
through  a systematic  study  by  a variation  of  the  following  parameters:  Debye 

numbers  of  0.01  to  1000,  speed  ratios  of  0.1  to  5.0,  non-dimensional  potentials 
of  -10  to  -200,  and  Knudsen  numbers  of  0.007  to  1000.  There  have  been,  how- 
ever, a number  of  assumptions  that  have  been  relaxed  or  modified  in  the  present 
study  which  will  be  discussed  in  detail.  For  example, 

0 the  effects  of  the  variation  of  the  electric  potential  on  the 
vehicle  have  been  investigated 

0 the  effects  of  vehicle  angle  of  attack  ■ 

0 the  consideration  of  internal  degrees  of  freedom 
In  addition  further  extensions  of  the  Monte  Carlo  method  including 

0 formulation  of  a procedure  to  describe  chemical  reactions 

0 preliminary  calculations  to  describe  the  collection  of 
negative  ions 


will  be  presented. 


» 

2.  POSITIVE  ION  COLLECTION  - A SYSTEMATIC  PARAMETER  STUDY 


The  vehicle  geometry  and  flow  configuration  shown  in  Figure  1 was  chosen 
to  represent  the  shapes  currently  in  use  to  house  rocket-borne  mass  spectrom- 
eters. The  front  face  of  the  cylinder,  the  orifice  plate,  was  assumed  to  be 
held  at  a negative  potential  with  respect  to  the  ambient  gas,  and  the  positive 
ion  flux  induced  by  the  resulting  electric  field,  as  well  as  the  density, 
temperature,  and  velocity  of  the  ions  in  the  flowfield  about  the  cylinder  were 
computed.  The  ambient  gas  was  assumed  to  be  composed  of  neutral  molecules, 
positive  ions  and  electrons.  The  concentration  of  ions  and  electrons  relative 
to  neutrals  is  actually  so  low  in  the  ionosphere  that  the  former  may  be  con- 
sidered trace  species.  Thus,  the  flow  of  the  charged  species  will  be  affected 
by  the  neutrals,  but  not  vice  versa.  Likewise,  binary  encounters  between 
charged  particles  can  be  neglected  and  the  only  charged  particle  interaction 
arises  through  the  electric  field. 

The  quantity  of  primary  concern  is  the  flux  of  positive  ions  to  the 
stagnation  point  of  the  cylinder  orifice  plate.  The  net  current  is  not  of 
interest;  hence  the  electrons  in  the  gas  are  o^ interest  insofar  as  their 
spatial  distribution  affects  the  net  space  charge.  The  electrons  were  pre- 
scribed to  be  an  inviscid  gas  in  Boltzmann  equilibrium  under  the  electric 
field.  This  assumption  was  shown  to  be  accurate  by  Vogenitz  (1973)  for  the 
positive  ion  collection. 

The  body  geometry  considered  in  this  study  was  a right  circular  cylinder 
at  zero  angle  of  attack.  The  electric  potential  of  the  face  of  the  cylinder 
with  respect  to  the  ambient  gas  was  specified  as  an  input  parameter,  whereas 
the  side  of  the  cylinder  was  assumed  to  be  at  zero  potential  with  respect 
to  the  ambient  gas.  This  means  that  the  small  charge  buildup  at  zero  current 
on  the  actual  flight  article  (the  floating  potential)  was  assumed  to  be  zero. 
Only  that  portion  of  the  cylinder  side  was  simulated  which  was  required  to 
compute  accurately  ion  current  to  the  front  face. 

The  hard-sphere  interaction  law  was  used  for  neutral-neutral  and  neutral- 
charged  particle  interactions.  The  numerical  magnitudes  of  collision  cross 
sections  for  neutrals  and  ions  were  taken  to  be  identical  in  the  present  study 
and  the  ion  and  electron  temperatures  in  the  ambient  gas  at  infinity  were 
assumed  equal  to  the  neutral  temperature.  The  method,  however,  could  treat 
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arbitrary  values  of  these  parameters  which  are  input  variables.  The  effects 
of  varying  the  speed  ratio,  front  face  potential  and  the  effects  of  different 
cross  sections  were  discussed  by  Sugimura  and  Vogenitz  (1975). 

In  Figure  2,  the  range  of  Debye  numbers  shown  as  a function  of  Knudsen 
number  or  altitude  was  taken  from  Marcisi  (1970).  This  figure  shows  the 
expected  Debye  number  range  for  a given  altitude  for  a vehicle  with  a 15  cm 
diameter.  As  mentioned  earlier  a systematic  variation  of  parameters  has  been 
completed  in  the  present  study  and  predessor  contracts  for  the  following  values 
and  has  been  reported  by  Sugimura  and  Vogenitz  (1975)  and  Sugimura  (1976). 

A 

0.007  1 -J  1 1000 

0.01  1000 

0.1  _ S < 5.0 

-200  < i < -10 
0 

The  Knudsen  number,  = 0.007,  corresponds  to  an  altitude  of  70  km,  the 
lowest  altitude  of  practical  interest.  % 

2.1  FREE  MOLECULE  FLOW  RESULTS 

The  ion  flux  coefficient  at  the  stagnation  is  shown  in  Figure  3 as  a 
function  of  Debye  number  for  speed  ratios  from  1.0  to  5.0.  It  can  be  seen 
that  Cj  is  constant  as  the  Debye  number  is  decreased  from  1000  to  a value 
between  0.10  and  1.0.  This  indicates  that  both  values  of  the  Debye  number 
lead  to  a Laplace  solution  of  the  Poisson  equation.  As  the  Debye  number  is 
decreased  further,  a power  law  relationship  between  Cj  and  \q/D  is  established. 

? The  variation  of  Cj  can  be  approximated  by  two  straight  lines  on  the  semi-log 

graph.  One  with  zero  slope  for  large  Debye  numbers  and  the  second  with  a 
slope  which  depends  on  the  speed  ratio.  This  simplification  will  overestimate 
slightly  at  the  transition  point  where  the  slope  of  the  curve  changes  but 
it  provides  a method  for  reducing  the  number  of  calculations  required  to  determine 
the  Debye  number  variation.  Also  shown  in  the  figure  is  the  flux  coefficient  if 
the  front  face  potential  is  taken  to  be  zero.  Since  the  flux  coefficient  Cj  is 
normalized  by  the  "freestream"  flux  this  value  is  unity  for  the  free  molecule 
flow  limit  but  is  less  than  unity  in  the  transition  regime. 

I' 
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STAGNATION  POINT  FLUX  COEFFICiLNi  AS  A FUNCTION  OF  DEBYE  NUMBER  FOR  VARIOUS  SPEED  RATIO 


FIGURE  3 


I 
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In  Figure  4 the  variation  of  Cj  with  speed  ratio  is  shown  for  Debye 
numbers  from  0.01  to  1000.  Again  very  little  difference  in  Cj  is  found  for 
Debye  numbers  1000  and  1.0.  For  speed  ratios  greater  than  unity  a power  law 
variation  exists  between  Cj  and  S with  the  exponent  depending  on  the  Debye 
number,  ^p/D.  The  values  of  Cj  increase  rapidly  for  speed  ratios  less  than 
unity,  due  mainly  to  the  fact  that  Cj  is  normalized  by  the  freestream 
velocity. 

The  previous  calculations  were  all  performed  with  a value  of  the  front 
face  potential  equal  to  -200.  In  Figure  5 the  stagnation  point  ion  flux 
coefficient  is  shown  to  increase  linearly  as  the  front  face  potential  is 
varied  from  .^  = 0 to  4^  = -200  for  a speed  ratio  ot  1.5.  These  detailed 
■-esults  for  a Tebye  number  of  0.1  v/ere  presented  previously  (Sugimura  and 
. '^enitz  1973;  If  a linear  variation  of  Cj  and  is  assumed  for  the  other 
•'bj. e numbers;  0.01,  1.0,  and  1000,  the  dashed  lines  in  Figure  5 result. 

The  'on  i.  • variation  across  the  front  face  of  the  collectinc  disk 
shown  in  : ' n'e  C fer  a Knudsen  number  of  ' /D  = 10  for  Oebye  numbers 
I G.IO  t’  Also  shown  in  cm  s fioure  i;  ■ result  for  a Kradsen 

' anH  . nir-bec  both  equal  ■;  TCO.  It  c:.n  be  seen  that  the  two 
SI  ''.s,  '•  , -GJO  - id  10,  for  u bebye  nuh'bi-r  of  1000  are  very  similar  and 
ese'.t  : ; ..lS  ’ 1 1 s ion  1 ess  and  nearly  ;:ol  1 isionless  cases,  f-'  ■ flux  for  ' 

/ ■ = !•  is  s'inhtly  lower  tnan  the  flux  for  /D  = lOOC,  which  si)ows  the  ; 

•.'A  of  a fir, tie  . .,mber  of  collisions.  However,  the  decrease  in  the  flux  | 

fficient  at  value  of  r^/R  = 0.07  is  due  to  a small  sample  size,  at  least  |i 

a factor  ten  less  than  the  sample  size  over  the  entire  front  face,  and  i 

the  stagnation  point  flux  is  found  by  extrapolating  through  the  other  points.  h 

A typical  sample  size  at  the  stagnation  point  was  approximately  3000  particles.  ;| 

i ; 

The  results  for  a Debye  number  of  1.0  and  1000  show  little  difference,  indi-  i| 

eating  again  that  the  Debye  number  range  between  1.0  and  1000  results  in  a !j 

Laplace  solution  of  Poisson's  equation.  As  the  Debye  number  is  reduced  fur-  j 

ther,  the  ion  flux  coefficient  also  decreases.  For  all  of  these  results,  | 

the  linear  variation  of  the  ion  flux  with  radius  which  characterizes  the 
near-free  molecule  flow  is  found  again.  | 

In  summary,  the  variation  of  all  the  important  parameters  for  the  collision- 
less case  are  represented  on  Figures  3 through  6.  The  speed  ratios  from  0.1 
through  5.0  and  Debye  numbers  from  0.01  through  1000  appear  to  cover  the  values 
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of  practical  concern.  The  range  of  the  potential  from  zero  to  -200  is  probably 

insufficient;  but,  if  the  linear  relationship  does  exist  between  Cj  and 

for  all  other  parameters  fixed,  the  parametric  study  for  the  collisionless 

limit  is  complete.  This  information  can  now  be  used  to  predict  the  ion  flux 

coefficient  Cj  for  any  combination  of  values  of  Debye  number,  ^q/D,  speed 

ratio,  S^.,  and  front  face  potential,  in  free  molecule  flow.  It  will  be 

shown  later  that  the  free  molecule  flow  result  is  valid  for  a Knudsen  number 

as  small  as  ^/D  = 1 which  corresponds  to  an  altitude  of  approximately  100  km.  | 

2.2  NEAR  FREE  MOLECULE  FLOW  | 

The  ion  flux  distribution  with  radius  for  a Knudsen  number  of  1.0  is 
shown  in  Figure  7.  By  comparing  these  results  with  those  for  a Knudsen  | 

number  of  10,  it  can  be  seen  that  the  effect  of  collisions  has  reduced  the 
ion  flux  by  approximately  one-third.  Again  for  large  Debye  numbers  the  flux 
variation  across  the  plate  is  about  a factor  of  2,  vihereas  for  the  small  Debye 
numbers,  e.g.,  ^^/D  = 0.05,  the  flux  is  approximately  constant  with  radius. 

This  indicates  that  the  sheath  is  very  close  to  the  body  surface  for  small 
Debye  numbers  and  the  electric  field  is  approximately  constant. 

2.3  TRANSITION  FLOW  (90  KM) 

The  density  variation  for  both  ions  and  neutrals  for  a Knudsen  number 
of  0.20  is  shown  in  Figure  3 for  Debye  numbers  0.02  and  1000.  It  is  seen 
that  there  is  a considerable  separation  between  the  ions  and  neutrals  in 
contrast  to  the  result  at  70  km  discussed  later.  There  is  also  no  indication 
of  a gasdynamic  shock  in  the  neutral  particles  which  is  expected  in  the  con- 
tinuum case. 

Figure  9 shows  the  flux  variation  across  the  front  face,  where  again 
the  stagnation  point  flux  is  determined  by  extrapolation.  In  extrapolating 
by  eye,  the  values  near  the  stagnation  have  less  weight  because  of  the  smaller 
number  of  samples  at  these  interior  points.  Similar  to  other  Knudsen  numbers, 
the  flux  becomes  uniform  as  the  Debye  number  decreases. 

The  variation  of  the  electric  potential  along  the  stagnation  line  is 
shown  in  Figure  10.  It  is  seen  tnat  for  Debye  numbers  greater  than  0.5  the 
DOtential  distrioution  does  not  change  and  is  given  by  the  Laplace  solution. 

For  smalle>"  Debye  numbers,  the  potential  decreases  faster  and  the  sheath 

becomes  smialler.  For  a Deoye  numoer  o^'  0.02,  the  potentia’'  becomes  positive  1 

at  approximately  7/-  - 2.5. 
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2.4  NEAR  CONTINUUM  FLOW  (70  m) 


Figure  11  shows  the  density  distributions  along  the  stagnation  line  for 
both  ions  and  neutral  particles  for  two  extreme  values  of  the  Debye  numbers, 

0.03  and  1000,  for  a Knudsen  number  of  x/D  = 0.007.  It  should  be  emphasized 
that  the  abscissa  is  measured  in  mean  free  path  units  X/A^  so  that  a value 
of  X/\^  = 20  is  in  terms  of  the  diameter 

X/D  = (X/A^)(-^/D)  =(20)(0.007)  = 0.14 

It  is  seen  from  Figure  11  that  the  first  effects  of  the  body  on  the  neutral 
gas  particles  occurs  at  approximately  X/a^  - 20  or  X - 0.14D. 

In  Figure  12  the  potential  field  variation  along  the  stagnation  line 
are  shown  for  the  two  Debye  numbers.  It  is  seen  that  in  both  cases  the 
electric  field  penetrates  far  upstream  to  accelerate  the  positive  ions  toward 
the  collecting  plate.  The  ions  are,  thus,  affected  far  upstream  of  the  neutral 
particle  gasdynamic  shock  identified  by  the  rise  in  density  from  X/-^  - 20  to 
X/\  - 10.  In  Figure  12,  large  differences  in  electric  potential  are  shown 

for  Debye  numbers  of  0.03  and  1000.  The  linear  behavior  of  the  potential  for 
■j^/D  = 1000  is  because  the  upstream  boundary  is  located  at  approximately 
X/  ^ - 60  which  corresponds  to  X - 0.42D.  For  this  case,  the  solution  of  the 
Poisson  equation  is  the  Laplace  solution  along  the  stagnation  line  determined 
for  a conf iguration  similar  to  the  field  between  two  parallel  plates  which 
gives  a linear  variation.  Although  the  ion  density  is  higher  near  the  body 
for  the  larger  Debye  number  the  electric  field  is  lower  and  the  net  effect 
on  the  ion  flux  at  the  stagnation  point  is  that  it  remains  approximately  con- 
stant, Cj  - 2.7,  for  the  Debye  number  range;  0.03  - Aj./D  - 1000. 

In  Figure  13  the  ion  flux  variation  across  the  front  face  shows  two 
effects.  First,  the  variation  with  Debye  number  in  the  range  0.03  - >|/D  - 1000 
is  very  small,  and  the  flux  at  the  stagnation  is  approximately  constant.  This 
result  is  shown  later  in  Figure  15.  Secondly,  the  ion  flux  peaks  sharply  at 
the  stagnation  point  reflecting  the  very  high  densities  for  this  near-continuum 
flow  field. 

Figure  14  shows  velocity  field  along  the  stagnation  line  for  both  ions 
and  neutrals.  Again,  as  with  the  density,  the  effect  of  the  electric  field 
upstream  of  the  disturbance  in  tne  neutral  molecules  is  apparent,  the  field 
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ION  VELOCITY  ALONG  STAGNATION  LINE 


FIGURE  14 


corresponding  to  the  larger  Debye  number  again  penetrating  further  from  the 
body.  It  should  be  pointed  out  that  the  ion  velocity  at  the  body  surface 
(X/\^  = 0)  is  not  zero,  unlike  the  velocity  for  the  neutrals,  since  the  ions 
are  assumed  to  be  converted  to  neutrals  when  they  strike  the  body.  These 
calculations  at  70  km  required  approximately  25  minutes  on  the  CDC  6600 
computer  for  the  most  time  consuming  calculations  at  small  Debye  number.  The 
ion  flux  samples  for  a typical  case  at  small  Debye  number  were  about  400  for 
the  stagnation  point  and  4000  over  the  entire  front  face. 

2.5  SUMMARY 

A surmary  of  the  results  for  the  parametric  study  is  shown  in  Figure  15, 
where  the  stagnation  ion  flux  coefficient  is  shown  as  a function  of  Debye 
number  for  Knudsen  numbers  equal  to  \,/D  = 0.007,  0.20,  1.0,  10.,  1000,  a 
speed  ratio  of  3.10,  and  a nondimensional  front  face  potential  of  -200.  Also 
included  on  this  figure  on  the  left  hand  side  are  the  computed  values  of  the 
flux  coefficient  in  the  zero-field  limit  ( i = 0),  which  represent  the  minimum 
value  for  a given  Knudsen  number.  Note  that  the  value  of  the  zero-field  limit 
for  the  free  molecule  case  (',/D  = ®)  is  unity  as  required. 

It  can  be  seen  that  the  variation  of  the  flux  coefficient  Cj  with  Debye 
number  exhibits  a common  trend  for  Knudsen  numbers  -0.20.  That  is  the  flux 
coefficient  is  constant  as  the  Debye  number  is  reduced  until  a "break-point" 
is  reached,  after  which  the  flux  coefficient  decreases  with  decreasing  Debye 
number.  A further  observation  is  that  the  break-point  occurs  at  a Debye 
number  approximately  equal  to  the  Knudsen  number  for  transitional  flows; 
e.g.  JD  - 0.20.  For  the  collisionless  flow  a /D  - 10,  the  break  point  occurs 
at  a value  of  >g/D  = 0.5  or  v.  body  radius.  The  effects  of  collisions  can 
oe  seen  in  the  results  for  Knudsen  number  unity,  e.g.  the  reduction  in  the 
^lux  coefficient  below  the  collisionless  limit.  However,  the  "break-point" 
remains  at  approximately  0.5.  The  results  for  the  near-continuum  flow  with 
Knudsen  number  0.007  exhibit  no  "break-point"  for  the  range  of  Debye  numbers 
considered.  If  the  "break-point"  is  assumed  to  occur  at  a Debye  number  approxi- 
mately equal  to  the  Knudsen  number  the  transition  should  occur  at  a Debye 
number  near  0.007.  Calculations  for  Debye  numbers  less  than  0.03  for  these 
conditions  have  net  been  successful  to  date  because  of  numerical  instabilities. 
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Although  these  results  do  not  represent  a complete  parametric  study  the 
following  trends  seem  reasonable. 

1)  The  flux  coefficient  decreases  with  decreasing  Debye  number  for 
a fixed  Knudsen  number. 

2)  The  flux  coefficient  decreases  with  decreasing  Knudsen  numbers  for 
fixed  Debye  numbers  greater  than  the  body  radius,  i.e.  a^^/D  ^ 0.5. 

3)  The  "break  point"  for  the  variation  of  the  ion  flux  coefficient 
with  Debye  number  can  be  estimated  by 

(Ad)bp  ^ minimum  R) 

i.e.  the  "break  point"  is  given  by  or  (^q)bp  = R which- 
ever is  less.  Therefore,  for  ^ determined 

by  taking  ^ 

These  limits  for  the  "break  point"  are  shown  on  Figure  2.  It  can 
be  seen  that  in  the  "range  of  interest"  below  90  km  the  ion  flux 
can  be  determined  by  the  infinite  Debye  number  limit. 

In  many  applications  the  quality  of  the  flux  entering  the  instrument  is  just 
as  important  as  the  number  of  particles.  Figure  16  gives  the  percentage  of 
free  stream  ions  in  the  stagnation  point  flux  as  a function  of  Knudsen  number. 
These  results  are  approximately  independent  of  Debye  number  for  the  range 
0.01  - q/D  - 1000.  This  shows,  for  example,  that  below  an  altitude  of  90  km 
( >^/D  = 0.2)  the  entire  flux  has  been  affected  by  the  body  and  only  for 
Knudsen  numbers  greater  than  10  is  the  flux  completely  undisturbed. 

Figure  17,  which  is  reproduced  from  Sugimura  and  Vogenitz  (1975)  shows 
the  variation  of  the  stagnation  point  flux  for  the  zero  field  limit  ( = 0) 

as  a function  of  Knudsen  number  or  altitude.  At  an  altitude  of  70  km 
= 0.007)  this  ion  flux  coefficient  is  approximately 


Cj  - 0.17 

Sonin  (1967)  gives  a continuum  prediction  for  the  zero  field  limit 


Cj  - 2 — 


-yil  S^Rg 
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Re  = Reynolds  number  = 


pQ  = Stagnation  density 
p = Free  Stream  density 

The  diffusion  coefficient  can  be  estimated  by  the  simole  relation  given  by 
Jeans  (1940). 

If  the  stagnation  conditions  are  taken  to  be  the  conditions  behind  a normal 
shock  for  a gas  with  a ratio  of  specific  heats  equal  to  1.4  the  value  of 
the  flux  coefficient  becomes 

Cj  = 0.16 

which  compares  favorably  with  the  Monte  Carlo  results.  However,  if  the 
velocity  and  temperature  are  fixed  the  variation  of  the  flux  coefficient 
with  mean  free  path  from  Sonin's  result  becomes 

Cj  'V  /T 

This  variation  shown  in  Figure  17  appears  reasonable  for  lower  altitudes 
but  cannot  be  extended  to  much  higher  altitudes.  A more  detailed  inves- 
tigation is  required  to  determine  the  correlation  between  the  continum 
prediction  and  the  Monte  Carlo  results.  Also  shown  on  Figure  17  is  the 
quality  of  the  flux  for  the  zero  field  case.  Here  the  flux  at  90km  is 
found  to  be  composed  entirely  of  ions  which  have  suffered  collisions 
with  neutral  molecules.  It  should  be  noted  here  that  the  collision  cross 
sections  for  ion-neutral  and  neutral-neutral  encounters  were  taken  to  be 
equal.  The  effects  of  different  cross  sections  are  discussed  by  Sugimura 
and  Vogenitz  (1975). 


STAGNATION  POINT  ION  FLUX  VARIATION  WITH  KNUUSEN  NUMBER  FOR  ZERO  PLATE  POTENTIAL 


FIGURE  17 


From  the  results  of  the  parametric  study,  a functional  relation  of 
the  form 


C]  " f(S..  >./D. 

can  be  determined  for  the  altitudes  which  are  characterized  by  Knudsen 
numbers  in  the  transitional  regime  or  larger,  > JD  = 0.20  to  in  the 
present  stuoy.  For  Knudsen  numbers  less  than  transitional  the  results  shown 
in  Figure  15  indicate  that  the  correlational  function  will  be  complicated 
by  the  location  of  tne  "break  point."  This  problem  becomes  unimportant, 
however,  if  the  range  of  parameters  is  restricted  to  the  "range  of  interest" 
shown  in  Figure  2. 

Furthermore,  in  the  collection  of  large  positive  (or  negative)  ions 
the  number  of  collisions  that  the  ions  experience  before  entering  the  mass 
spectrometer  is  important  because  of  perturbations  to  the  ambient  states, 
i.e.,  chemical,  internal,  etc.  The  result  indicated  here  that  the  fraction 
of  free  stream  ions  is  independent  of  the  Debye  number  and  electric  field  may 
be  of  some  use  to  instrument  designers. 
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3.  SIDEWALL  POTENTIAL  EFFECTS 


ft 

1. 


In  the  Monte  Carlo  results  reported  previously  the  space  charge  effects  on 
the  electric  potential  of  the  sidewall  were  neglected.  That  is  the  boundary  con- 
dition for  the  potential  on  the  vehicle  side  walls  was  assumed  to  be  zero.  How- 
ever, it  is  known  that  typical  values  of  floating  potentials  are  1 to  3 volts 
negative.  In  order  to  investigate  the  effects  of  finite  space  charge  on  the 
sidewall  the  boundary  conditions  for  the  potential  were  modified  in  the  following 
way. 


The  applied  front  face  potential,  in  a non-dimensional  form  is: 


*0  “ tr  " T 


$„(volts) 


The  non-dimensional  sidewall  potential  is: 


- 11  ,600  j 


SW 


and  the  potential  at  the  corner  of  the  vehicle  is  defined  as: 

'c  i g (.‘o  + 

Recall  that  in  the  introduction  the  sysc(.:"atic  parameter  study  was  to  be  made  for 
a fixed  non-dimensional  potential  of 

= -200 

Therefore,  in  the  absence  of  a sidewall  potential  the  ambient  temperature  does 
not  enter  into  the  analysis.  When  the  space  charge  effect  on  the  vehicle  is 
included  by  assuming  a fixed  sidewall  voltage,  e.g.,  = -1  volt,  the  ambient 


temperature  becomes  a parameter.  If  a value  of  i: 
following  results  are  obtained: 


SW 


-1  volt  is  assumed  the 
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^ -IK  f -p<»*  ■ 


T 


A - n fi(!h  ’ 

Altitude  (km)  T(»K)  *SW''’®™T7’ 


70 

216 

-53.7 

90 

184 

-63.0 

no 

257 

•45.2 

130 

534 

-21.7 

Th<s  table  shows  that  the  minimum  temperature  occurs  at  90  km  and  gives  the 
maximum  sidewall  potential. 


One  of  the  reasons  for  examining  the  effects  of  the  sidewaii  poteniial 
buildup  was  that  the  solution  of  Poisson's  equation  for  the  electric  potential 
became  unstable  for  small  Debye  lengths.  This  difficulty  was  caused  by 
large  changes  in  the  electric  potential  over  distances  comparable  to 


r 


The  ion  flux  at  90  km  was  found  to  be  greater  than  at  70  km.  If  a sidewall  volt- 
age of  -1  volt  is  assumed  the  ion  ■^lux  at  90  km  is  less  than  at  70  km  because 
the  non-dimensional  sidewall  potential  is  larger  at  90  km.  Since  the  objectives 
of  this  phase  of  the  study  were  to  obtain  general  trends  in  the  variation  of 
the  parameters  it  was  decided  to  minimize  the  number  of  variables  by  ignoring 
the  sidewall  potential  in  most  of  our  calculations.  In  the  future  if  the  positive 
ion  flux  is  to  be  predicted  along  a given  flight  trajectory  the  front  face  and 
sidewall  voltages  should  be  fixed  and  all  the  non-dimensional  potentials 
and  should  be  a function  of  altitude. 
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4.  ANGLE  OF  ATTACK  EFFECTS 


The  difficulty  in  considering  angle  of  attack  effects  is  that  the  flov; 
symmetry  is  destroyed  and  a three  dimensional  flowfield  must  be  considered. 

A Monte  Carlo  simulation  code  to  simulate  the  flow  of  neutral  molecules  about 
general  body  shapes  at  arbitrary  angles  of  attack  is  currently  operational 
the  AFGL  computing  system.  Modifications  to  this  program  are  necessary 
before  any  charged  particle  calculations  can  be  made.  In  particular  a three- 
dimensional  Poisson  equation  solution  is  required  to  calculate  an  electric 
field  consistent  with  the  charged  particle  flowfield. 

The  geometry  for  the  body  at  angle  of  attack  is  shown  in  Figure  18.  The 
parameters  which  characterize  the  three  dimensional  effects  are  the  angle  of 
attack,  i,  and  the  azimuthal  angle,  o.  In  the  two  dimensional  case  a = 0 
and  there  are  no  variations  with  e. 

Figure  19  shows  the  density  distribution  of  neutral  molecules  in  a plane 
parallel  to  the  front  face  as  a function  of  the  azimuthal  angle  for  three  values 
of  the  angle  of  attack:  at  = 0,  10°,  and  20°.  This  distribution  is  given  for  a. 

ooint  located  in  a plane  a distance  0.091R  downstream  of  the  front  face  (R  is 
the  bcdy  radius  shown  in  Figure  18)and  a distance  0.62R  from  the  vehicle  sidewall 
As  expected  the  density  for  non  zero  angles  of  attack  is  hignest  for  small  e where 
tr^e  flow  is  compressed  and  lowest  for  6 = 180°  where  the  flow  has  expanded.  The 

density  changes  are  also  greater  for  higher  angles  of  attack.  An  interesting 
coincidence  in  this  case  is  that  the  density  for  all  three  angles  of  attack  are 
approximately  the  same  at  ■ = 100°.  The  value  of  the  dencity  ratio  at  the  stag- 
nation point  does  not  differ  by  more  than  10%  as  the  angle  of  attack  varies  from 
0°  to  20°.  Therefore,  the  flux  entering  the  orifice  at  the  stagnation  point  is 
only  slightly  affected  for  these  flight  conditions. 


COORDINATES  FOR  THREE  DIMENSIONAL  FLOWFIELD 


FIGURE  18 


AZIMUTHAL  VARIATION  OF  DENSITY 


FIGURE  19 


TEMPERATURE  DISTRIBUTIONS  ALONG  STAGNATION  LINE 


5.  INTERNAL  DEGREES  OF  FREEDOM 


All  of  the  results  presented  were  based  on  a hard  sphere  molecular  model 
with  no  internal  degrees  of  freedom.  A more  realistic  molecular  model  would 
include  the  effects  of  rotational  and  vibrational  degrees  of  freedom.  For 
example.  Figure  20  shows  a comparison  between  the  average  kinetic  or  transla- 
tional temperature  and  the  rotational  temperature  along  the  stagnation  line. 
This  calculation  was  performed  for  a neutral  gas  at  a Knudsen  number  of  0.10 
corresponding  to  an  altitude  of  88  Kiti  and  a speed  ratio  of  3.‘'0.  If  the 
gas  were  in  thermodynamic  equilibrium  the  two  temperatures  would  be  identical. 
However,  the  kinetic  temperature  is  shown  to  rise  very  rapidly  while  the 
rotational  temperature  remains  much  lower.  Finally  at  a distance  of 
approximately 


the  temperatures  are  again  in  eouilibnum.  A simple  explanation  of  this 
effect  can  be  made  by  notinc  t‘‘it  rotational  equilibrium  requires  many  more 
collisions  than  ^ ranslational  equilibrium,  e.g . typical ly  a ratio  of  10:1. 

The  flow  is  assumed  to  be  ir  eouilibrium  far  upstream  of  the  body  but  as 
the  presence  of  the  body  is  encounterec  the  density  increases.  However, 
this  density  rise  is  not  rufficient  to  maintain  equilibrium  and  the  *ota- 
tional  temperature  falls  below  the  kinetic  temperature.  Finally,  as  the 
density  increases  through  the  "shock"  a level  is  reached  where  the  collisions 
are  sufficient  to  reestablish  equilibrium.  Also  shov/n  in  Figure  20  is  the 
kinetic  temperature  distribution  if  the  rotational  energy  is  igiiored.  It 
is  seen  that  without  the  rotational  energy  which  acts  as  an  enemy  sink  the 
translational  temperature  is  uniformly  higher. 

A more  explicit  example  of  the  effect  of  molecular  collisions  i'.  giver 
in  Figure  21  where  the  stagnation  line  temperature  distribution  is  presented 
for  a Knudsen  number  of  0.007  corresponding  to  an  altitude  of  70  Km  and  . 
speed  ratio  of  3.10.  The  density  increase  and  the  collision  frequnm.ie  are 
iiuch  higher  than  for  the  results  at  a Knudsen  number  of  0.10.  The  effects 
of  tne  increased  collisioris  on  the  temperature  distributions  is  seer,  in  the 
near  equilibrium  in  the  two  temperatures  for  x/a^  -5.0.  I<ote  that  the  abscissa 
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FIGURE  21 


• 


in  both  Figures  20  and  21  is  the  distance,  x,  divided  by  the  mean  free  path, 

A^,  which  is  approximately  15  times  greater  at  88  Km. 

If  the  velocity  distribution  function  is  examined,  the  initial  equilibrium 
Boltzmann  distributuion  far  upstream  of  the  body  is  distorted  as  the  non 
equilibrium  effects  are  encountered.  Again,  as  the  collision  frequency 
increases  near  the  body  equilibrium  conditions  are  reestablished. 

The  effects  of  rotational  non  equilibrium  on  the  collected  ion  flux  are 
not  obvious.  However,  the  addition  of  internal  degrees  of  freedom  acts  as  an 
energy  sink  which  certainly  influences  the  other  flow  paran'eters.  From  the 
calculations  at  Knudsen  number  0.10  the  stagnation  point  flux  for  zero  potential 
= 0)  differs  by  about  10%  when  rotational  energy  is  included,  further 
computations  are  required  to  determine  the  extent  of  this  difference  due  to 
statistical  scatter. 
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CHEMICAL  REACTIONS 


It  has  been  shown  in  Section  2 that  the  properties  of  the  ions  collected 
by  the  mass  spectrometer  are  distorted  from  the  ambient  values.  In  particular 
for  large  highly  clustered  ions  the  distortion  during  sampling  may  be  serious 
enough  to  prevent  the  determination  of  the  true  cluster  ion  distribution. 

The  effects  of  the  change  in  chemical  composition  can  be  assessed  in  a straight- 
forward way  through  the  direct  simulation  Monte  Carlo  method  by  treating  the 
chemical  reactions  as  collision  processes. 

A typical  bimolecular  reaction  may  be  written  as 


A + B 


C + D 


where  A,  B,  C,  D represent  separate  species,  e.g.  molecules,  ions,  electrons. 
If  the  molecules  can  be  idealized  as  reactive  hard  spheres  with  total  elastic 
scattering  cross  section,  o^,  the  total  reaction  cro:s  section  is  given  by 


0 for  E < E* 


★ 

where  E = threshold  energy 

E = relative  energy  in  the  collision 
and  the  reaction  rate  coefficient,  k(T),  is  ■'^ound  to  be 


(T)  = 


1/2 

exp 


where  „ = reduced  niass. 

If  k(T)  and  E*  are  known  for  the  important  processes  in  the  ionosphere,  the 

cross  section  ri  can  be  determined  and  fed  directly  into  the  Monte  Carlo  Simulation 
0 

calculations.  This  model  can  be  extended  to  include  internal  degrees  of  freedom 
and  deviations  from  hard  sphere  interactions.  Although  no  calculations  involving 
chemical  reactions  has  been  attempted  for  the  ion  collection  processes  the  model 
has  been  used  successfully  to  predict  radiation  effects  in  high  altitude  plumes. 
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7.  NEGATIVE  ION  COLLECTION 


The  assumptions  and  the  complete  formulation  of  the  Monte  Carlo  simu- 
lation are  detailed  in  Vcgenitz  (1973).  Most  important  of  the  assumptions 
regards  the  treatment  of  the  electrons.  For  the  positive  ion  collection,  the 
front  face  is  negatively  charged  which  repels  the  electrons  and  attracts  the 
positive  ions.  The  negative  front  face  potential  was  assumed  large  enough 
so  that  the  electrons  could  be  taken  to  be  in  equilibrium  with  the  local 
potential  and  determined  by  the  Boltzmann  distribution. 


= exp(e^,/kTp) 


where 


= electron  number  density 
= free  stream  value 
= electric  potential 
= electron  temperature 
= charoe  on  electron 


This  assumption  was  shown  to  be  accurate  and  was  then  used  throughout  the 
previous  studies.  The  great  simplification  of  this  assumption  was  that  the 
electron  trajectories  could  be  ignored  and  only  the  ion  motion  v;as  computed 
with  tne  only  coupling  between  the  ion  and  electron  distributions  occurring 
phrough  the  electric  potential  determined  from  Poisson  s equation 


■ (^)  - '•e> 


where 


ti  = = non-dimensional  electric  potential 


X = mean  free  path 
= Debye  length 
'*i 

N = positive  ion  number  density  ratio 

i N. 

on 

— Mp  , 

N = = electron  number  density  ratio 

e Ng 
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The  collection  of  negative  ions  requires  a positive  collecting  potential 
on  the  front  face  which,  however,  also  attracts  the  electrons  which  will  reduce 
the  collecting  potential  and  the  effectiveness  of  the  instrument.  The  sampling 
of  negative  ions  requires  a circular  double  disk  configuration;  Sherman  and  Parker 
(1970),  on  the  front  face  to  repel  the  electrons  while  still  collecting  negative 
ions.  The  geometry  and  the  potential  distribution  along  the  front  face  of  the 
instrument  payload  will  be  approximated  as  shown  in  Figure  22.  Although  the 
face  potential,  to,  would  be  greater  than  the  mask  potential,  $e,  the  mask  has 
greater  area.  The  radii  r^  and  r^  define  the  disks  corresponding  to  the  two 
potentials.  Therefore,  only  the  proper  ratio  of  to  <1^  will  allow  the  ions 
to  pass  through  the  collecting  orifice  for  a given  set  of  free  stream  and 
boundary  conditions.  The  electron  motion  in  this  case  cannot  be  ignored,  and 
the  equilibrium  assumption  for  the  electrons  is  no  longer  valid.  This  is  the 
major  complication  introduced  by  the  negative  ion  collection  which  was  not 
encountered  in  the  sampling  of  positive  ions. 

These  complications  can  be  seen  by  comparing  the  electron  and  ion  veloci- 
ties. The  thermal  velocities  for  the  electrons  and  ions  are  given  by 


where 


for  ions 


for  electrons 


^i*  ^e 

m. , m 
1 e 


ion  and  electron  temperatures 
ion  and  electron  masses 


‘I-  = positive  collection  potential  for  negative  ion 
c 

$ = negative  electron  mask  potential 


Figure  22.  Front  Face  Geometry 


The  ratio  of  these  speeds  are  given  by 


which  fcr  equal  temperatures,  Ti  = Te,  and  an  ion  molecular  weight  of  30 
becomes 

Ve 

^ = 235 

V. 

Since  the  vehicle  or  free  stream  velocity,  U , is  roughly  3 to  5 times  the 
ion  the*-mal  speed,  the  ions  are  characterized  by  the  free  stream  velocity 
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Ui  - U - 4 V. 

1 

whereas  the  typical  electron  velocity  is  the  thermal  velocity 


Up  - - 235  V.  - 60  U, 

® e 1 


Therefore,  for  a gi.en  time  increment,  the  electron  travels  a distance  60 
times  the  ion  displacement.  If  the  time  step  in  the  Monte  Carlo  calcula- 
tions is  given  by  :ti  for  the  ions,  the  result  requires  60  electron  time 
steps  to  equal  /ti.  In  other  v/ords,  each  ion  calculation  requires  60 
electron  calculations.  This  number  can  be  reduced  by  recalling  that  both 
the  electrons  and  ions  are  trace  species  and  that  encounters  between  charged 
particles  are  neglected  except  through  their  coupling  in  the  electric  poten- 
tial. The  only  direct  interactions  are  collisions  between  the  charged 
particles  and  the  neutral  molecules  which  can  be  characterized  by  the  mean 
■^^ee  paths 


/2  Q N 
en  n 


Oin  "n 


for  electrons/neutrals 


for  ions/neutrals 


where 


N = number  density  of  neutrals 
n '' 

Q , Q.  = collision  cross  sections 
^en’  ^in 


Assuming  hard  spheres,  the  collision  cross  sections  'I'or  two  di'f'ferent  parti- 
cles is  given  by 


a-j,  C2  = particle  diameters 


1*4*  %• 
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For  the  ions,  o,.  - o„  and  Qj„  - tt  a„^,  but  for  the  electrons  o„  « o„  and 

^ 2 ’ ^ ^ e n 

Qen  - TT  Op^/4  The  ratio  of  mean  free  paths  is  then  given  by 


X Q. 
en  ^in  . 


Therefore,  if  an  ion  travels  a distance  in  time  increment  Ati  an  electron 
disolacement  can  be  x^^  4 x^^.  This  reduces  the  number  of  electron  time 

steps  to  equal  Ati  from  60  to  60/4  =15.  In  rough  terms,  this  means  that 
ttie  negative  ion  collection  calculations  will  take  approximately  15  times 
as  long  as  the  equivalent  positive  ion  calculation  which  ignored  the  electron 
trajectories. 

A procedure  to  be  used  in  the  negative  ion  sampling  involving  a complex 
iteration  procedure  has  been  formulated  and  is  being  implemented.  However, 
preliminary  results  for  the  limiting  case  of  free  molecule  flow  and  large 
Debye  number  are  presented  along  with  calculations  for  electron  trajectories. 

7.1  FREE  MOLECULE  FLOW  AND  INFINITE  DEBYE  NUMBER 

In  the  limit  of  very  large  Knudsen  and  Debye  numbers  the  collection  of 
negative  ions  can  be  predicted  with  the  same  computer  code  used  for  the 
positive  ion  calculations.  For  large  Debye  numbers  the  electric  field  is 
not  affected  by  the  ion-electron  distribution  and  the  potential  is  given  by 
the  Laplace  solution.  This  case  is  not  insignificant  since  it  was  found  in 
Section  2 that  the  Laplace  solution  was  valid  for  Debye  numbers  as  small  as 
1.0.  Although  the  validity  of  the  Laplace  solution  must  be  determined  again 
for  the  negative  ion  processes  it  should  be  valid  for  a similar  range  of 
Debye  numbers.  The  geometrical  configuration  for  the  collection  of  negative 
ions  with  an  electron  mask  was  shown  in  Figure  22.  The  Laplace  solution  for 
this  double  disk  geometry  was  given  by  Sherman  and  Parker  (1970)  and  can  be 
characterized  by  the  following  parameters: 

ratio  of  disk  radii,  a = r /r 

0 e 

ratio  of  potentials,  • = 


1 
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Since  the  attracting  potential,  is  positive  and  the  el,ectron  mask  potential 
is  negative  it  is  possible  that  the  potential  along  the  stagnation  line  can 
become  negative  and  repel  1 the  negative  ions.  Sherman  and  Parker  shov;  that 
the  critical  value  of  t is  given  by 


1 - V ^ ■ “ 

and  is  shown  in  Figure  23.  If  1y|  > |y^|  there  is  no  repelling  potential  and 

IyI  < IY(.1  the  potential  becomes  negative  somewhere  along  the  stagnation  line. 
The  variation  of  the  potential  along  the  stagnation  line  from  the  analytic  result 
will  be  compared  to  the  finite  difference  solution  used  in  the  Monte  Carlo 
iteration  procedure. 

The  Monte  Carlo  calculations  for  the  negative  ion  collection  in  the  collision- 
less and  large  Debye  number  limit  were  performed  for  the  following  choice  of 
parameters. 


~ = 1000 


= 1000 


S = 3.10 


a = ~ = 0.50 
' e 

f = 200 
c 


= -40  and  -20 


The  two  values  of  gives  values  of  y equal  to 


Y = -5  and  -10 


which  from  Figure  23  indicates  a repelling  potential  for  the  first  value  but 
not  for  the  second.  In  Figures  24  and  25  the  finite  di-^ference  results  along 
two  lines  r = 0 (stagnation)  and  r = V./ ^ are  shown  and  compared  to  the  analytic 
solution.  Figure  24  shows  excellent  agreement  for  the  two  results  along  the 
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0 

R/2 


ANALYTIC  SOLUTION 
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(SHERMAN  AND  4>  = 200  = R/2 

PARKER,  1970)  ° ° 


stagnation  line  where  both  show  a change  in  sign  of  the  potential  at  approxi- 
mately two  radii  downstream.  The  agreement  in  Figure  25  is  not  as  good  but 

this  discrepancy  is  probably  due  to  the  finite  node  spacing  in  the  numerical 

calculations  which  makes  the  definition  of  a somewhat  ambiguous.  The  results 
for  the  ion  flux  variation  across  the  front  face  in  Figure  26  show  the  reduction 

of  the  stagnation  point  flux  when  y is  changed  from  -10  to  -5.  If  the  value  of 

IyI  vjas  further  decreased,  a point  would  be  reached  where  the  negative  ion  flux 
at  the  stagnation  point  would  vanish.  Also  shown  on  the  same  figure  is  the  ^lux 
variation  of  positive  ions  with  the  same  parameters  except  for  a negative  attract- 
ing potential . 


= -200 

The  difference  in  the  negative  and  positive  ion  results  is  due  to  the  prescribed 
potential  on  the  front  face.  For  the  positive  ion  collection  the  potential  is 
negative  and  uniform  but  for  the  negative  ion  collection  the  electron  mask  is  placed 
over  the  mass  spectrometer  housing  to  retard  the  electron  current.  This  results 
in  a positively  charged  central  region  and  a negatively  charged  outer  disk. 

In  both  cases  the  ion  flux  is  not  the  only  contribution  to  the  measured 
current.  Electrons  may  enter  the  mass  spectrometer  in  both  cases  but  will 
be  more  likely  to  do  so  in  the  collection  of  negative  ions. 

This  same  procedure  for  large  Debye  numbers  can  be  carried  out  for 
Knudsen  numbers  at  which  collisions  must  be  taken  into  account  and  represents 
one  0^  the  simpler  options  of  the  general  iteration  procedure.  A parametric 
study  for  this  limiting  case  will  provide  valuable  information  for  the  gen- 
eral process  of  negative  ion  collection  without  the  time  consuming  calcula- 
tions required  to  describe  the  electron  motion. 

7.2  ELECTRON  MOTION 

In  order  to  assess  the  computer  time  requirements  and  accuracy  of  Monte 
6-arTo  method  for  the  prediction  of  the  electron  motion,  the  following  calcu- 
lations were  made.  Figure  27  shows  the  electron  distribution  along  the  stag- 
nation line  for  the  following  conditions; 
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Collisionless  flow 


Large  Debye  number  >>  1 

Negative  Potential  4.  = -200  (uniform  across 

° front  face) 

In  the  calculations  for  the  collection  of  positive  ions,  the  electrons  were 
assumed  to  be  in  Boltzmann  equilibrium  for  large  negative  front  face  potentials. 
For  the  calculations  shown  in  this  figure  the  Monte  Carlo  calculations  for  the 
electrons  are  compare  favorably^ with  the  results  of  the  equilibrium  electron 
gas  assumption.  The  electron  speed  ratios  for  the  calculations  were  taken 
to  be  0.10  and  0.014  corresponding  to  ion  speed  ratios  of  24  and  3.4, 
respectively.  The  two  results  for  these  speed  ratios  are  similar  enough  to 
conclude  that  for  very  small  speed  ratios  the  electron  motion  is  independent 
of  speed  ratio  for  a repulsive  potential.  The  increased  scatter  in  the 
results  for  a speed  ratio  of  0.014  is  due  to  the  smaller  movement  time 
interval  for  a fixed  time  step,  i.e.,  the  calculation  for  the  smaller  speed 
ratio  represents  a shorter  simulation  time.  Since  the  larger  speed  ratio 
calculation  converges  faster,  it  would  be  useful  to  determine  the  largest 
speed  ratio  for  the  given  conditions  which  results  in  a solution  independent 
of  the  speed  ratio. 

Another  series  of  Monte  Carlo  calculations  were  made  to  determine  the 
requirements  to  predict  the  electron  motion  for  a positive,  attractive  front 
face  potential.  For  these  calculations  the  flow  was  again  assumed  collision- 
less, the  Debye  number  large  enough  for  the  Laplace  solution  to  be  valid  and 

the  front  face  at  a non  dimensional  potential  equal  to  200.  No  electron  mask 
was  placed  over  the  front  face  and  the  purpose  of  these  calculations  was  to 
properly  simulate  the  electron  motion  with  an  attractive  field.  Figures  28 
through  30  show  the  potential  distribution,  number  density  and  velocity  along 
the  stagnation  line.  In  the  collisionless  case,  the  electrical  energy  and 
kinetic  energy  of  the  electron  are  directly  related  as  shown  in  Figures  28 
and  30.  A number  of  calculations  were  made  for  various  values  of.  Atm,  the 
movement  time  interval.  The  value  shown  on  the  figures  was  found  to  be  close 
to  optimum  between  the  larger  values  which  led  to  convective  errors  and  smaller 
values  which  required  excessive  computer  time  to  reach  a steady  state. 


51 


VARIATION  OF  ELECTRIC  POTENTIAL 


8.  CONCLUSIONS 


The  results  of  the  systematic  parameter  study  for  the  collection  of  positive 
ions  were  discussed  in  Section  2.  Most  of  the  important  collection  variables 
were  investigated  over  appropriate  ranges  of  interest.  Additional  parameters 
such  as  the  effects  of  vehicle  angle  of  attack  and  the  molecular  internal 
degrees  of  freedom  were  examined.  However,  only  a few  calculations  were  made 
to  assess  the  effects  of  relaxing  the  original  assumptions,  e.g.,  zero  angle 
of  attack  and  neglect  of  the  internal  degrees  of  freedom.  Results  of  the 
angle  of  attack  calculations  indicate  that  the  flux  coefficient  is  only 
slightly  affected,  on  the  order  of  10%,  for  moderate  angles  of  attack,  e.g. 
less  than  20  degrees.  The  effects  of  rotational  energy  was  also  formed  to  be 
about  a 10%  effect  for  the  flow  conditions  considered.  Further  work  for  the 
positive  ion  collection  should  include  a systematic  study  of  these  effects 
and  the  role  of  chemical  reactions  on  the  collection  process. 

The  formulation  for  collection  of  negative  ions  has  been  completed  and  the 
preliminary  calculations  indicate  that  the  electrons  are  described  properly. 
Calculations  completed  for  the  negative  ion  flux  in  the  collisionless  limit 
for  large  Debye  numbers  will  be  continued  and  extended  into  the  transition 
regime. 
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